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The paper describes the multiple rotating frame technique for designing modulated rf- 
fields, that perform broadband heteronuclear decoupling in solution NMR spectroscopy. The 
^ CT^ decoupling is understood by performing a sequence of coordinate transformations, each of 
^ which demodulates a component of the Rf-ficld to a static component, that progressively 

^ averages the chemical shift and dipolar interaction. We show that by increasing the number 

^ of modulations in the decoupling field, the ratio of dispersion in the chemical shift to the 

^ strength of the rf-field is successively reduced in progressive frames. The known decoupling 

O methods like continuous wave decoupling, TPPM etc, are special cases of this method and 
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their performance improves by adding additional modulations in the decoupling field. The 
technique is also expected to find use in designing decoupling pulse sequences in Solid State 
NMR spectroscopy and design of various excitation, inversion and mixing sequences. 
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1 Introduction 



Heteronuclear decoupling methods have a long history in NMR spectroscopy [H El |3l HJ 

EiEiEiHiiHiiTniiniiiaiHiEiEiiiniiniiTHiiiniE^ 

EO] and in vivo applications [211 E21 E31 EH ESI ES] • The goal of broadband heteronuclear 
decoupling sequences is to observe a spin S by irradiating spin I that is coupled to S in 
order to simplify the spectra and to increase the signal-to-noise ratio. At the same time, the 
decoupling sequence should introduce only a minimal amount of artifacts, such as decoupling 
sidebands. Furthermore, in order to avoid undesirable sample heating or damage to the 
probe, the radio frequency (rf) power of the decoupling sequence should be as small as 
possible. This is of particular importance in medical imaging or in vivo spectroscopy of 
humans [3ll E2l ESI El ES EH] . 

The earliest heteronuclear decoupling methods were based on cw irradiation and 
noise decoupling [12J. Significantly improved decoupling sequences were found based on 
composite [8l[9l[T3l[Tl[T6l[I71[l8l[l9]or shaped [211 [251 [261 [271 |28l EO] inversion pulses in 
combination with highly compensated cycles and supercycles [HI |9l [161 E3, [Ml [391 Hffl [ID [121 
Theoretical approaches that have been used for the analysis and design of decoupling 
sequences include average Hamiltonian [151 HH SS] and Floquet [23 EO] theory. 

Here, we introduce a powerful approach for the design of heteronuclear decoupling se- 
quences based on multiple rotating frame technique. 

In the absence of chemical shifts, the most straightforward way to decouple two heteronu- 
clear spins / and 5* is to just irradiate with full power on resonance to spin /. The presence of 
large chemical shifts makes this not the best strategy as using a CW irradiation, one creates 
an effective field that is not perpendicular to the coupling interaction and therefore part of 
the coupling interaction, parallel to this field is not averaged. We show that it possible to de- 
sign a multiply modulated rf-field, whose effect is best understood by performing a sequence 
of coordinate transformations, specific to each chemical shift. Each transformation reduces 
the chemical shift interaction and averages part of the coupling until the chemical shifts 
are made arbitrarily small. Each transformation demodulates a component of the multiply- 
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modulated rf-field to a static component. The subsequent transformation is in the rotating 
frame defined by new static component and the residual chemical shift. We show that by 
increasing the number of modulations in the decoupling field, we can significantly improve 
the decoupling. Methods to make this method robust to rf-inhomogeneity are discussed. 

2 Heteronuclear Decoupling 

Consider two heteronuclear spins / and S. The Hamiltonian of the spin system takes the 
form 



where oji and ojs are the chemical shifts of spin I and S. Assuming spin S is being 
observed, we can without loss of generality take the chemical shift of the spin S to be zero 
and denoting uj = cuq, consider the natural Hamiltonian 



Hq = oojlz + 2itJIzSz + oJsS; 



(1) 



Ho = (jJoI^+ 2jrJIzSz^, 

Hcs He 

where ujq G [— Co,Co]. Now, consider the rf-Hamiltonian 



(2) 



Hi 



= wqIx + {wi sin vit + W2 cos vit sin V2t + . . .)Iy. 



(3) 



A more general form of the rf-Hamiltonian is then 



N k-l 



h'q = + {^Wk n 



cos Vjt sin^jfct)/. 



(4) 



k=l j=l 



Ao{t) 



where Vk > Vk+i- We adopt the notation Wk — 2 



'wk, and rewrite 



Ho = + Ao{t)Iy + 27iJIzS, 



Z-i 



(5) 
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where 2:1(0;) is the unit vector along the direction uiz + wqx and define 9i{ijj) — tan~^(^), 
the spread of frequencies u e [wq, ^Jwq-\- Cq] . Note, zi{u), is different for each u. Now, by 
choosing vi in the first transformation as exactly the center of this spread and transforming 
into a frame rotating around /^i(aj), with frequency Di, we get the Hamiltonian 



//i = (a) - vi) 7,,(^) + 27rJcos^i(a;)4,(^)5, + wj^, + A^{t)Iy +H';{t) + iff, (6) 

^ -V^ ^ ^ V ^ 

where H[ is the demodulated part of the Rf-Hamiltonian iJg, in the interaction frame of 
Izi{w), that resembles Hq by design. H'{{t) and H'{'{t) are the fast oscillating parts of the rf 
and coupling Hamiltonian that we ideally want to average out and we neglect these terms for 
now. H"'{t) is simply the part of the coupling perpendicular to the effective field direction 
Izi{u>), that oscillates with frequency vi. The new frequency fi{uj) e [— ci, ci], where ci < cq. 
Ai, H'( and H'C are written in their general form below. The system obtained after first 
coordinate transformation has the desired feature that the ratio of chemical shift spread to 
Rf-strength ai = ^^is reduced over cto = ^ for the original system. For example, if cq = wq, 
then ^ = v^- 1. 

We can now iterate the above construction. Wc go into successive rotating frames around 
axis Izk{(xi) with frequency Vk- Unit vectors {xk{<jj),yk{oj), Zkioj)) define the k^^ frame, where we 
suppress the argument ui subsequently, /^(a;) is the chemical shift in the k*^ rotating frame, 
starting with ui in Hq. represents the limit of the chemical shifts in the /c*'* coordinate 
frame and Wk — '2.~^Wk is the strength of the rf- field along the direction Xk- 

Hk = + ^4. + 2nJkh,Sz + A,{t)Iy + H'l{t) + H'^{t), (7) 

n rri—l 

^feW = ^{ ^™ n COs{Vit)s\Ti{Vmt)}Iy, (8) 

m=k+l i=k+l 

Hl{t) = exp{z2vklz,t){-'^h, + Ak{t)Iy) exp{-z2vkh,t), (9) 
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H'l'it) = Jk exp{ivkIzkt)IxkSz ewi-i^kh^t) (10) 



^Cl + Wl-Wk ^Cl + Wl + Wk ^ Wk 

Ck+i = - ; Vk+i = - ^ ; tan6/fc(a;) = — — - (11) 

n 

Jfc= ncos^,.(a;)J. (12) 

k=\ 

As discussed in the next section the iterated relations on insures that the ratio ak = ^ 
is decreasing and 6*^ — )• f . This ensures that — 0. 



3 Scaling 



To fix ideas, we choose Wk = Wq, i.e., Wk = 2 ^wq. From equation (11), we have the relation, 



Uk+i = ^ I + al - I < Uk. (13) 
This ensures that ak is decreasing. We explore two limits, for <^ 1, 



ttfc 2 

afc 1, we have 



(14) 



ttfc+i ~ ttfc - 1. (15) 

Since — )■ 0, we have 6'^ — |. This ensures that — 0. 

We evaluate the root mean square amplitude for this rf-field as — )■ cxd 



Aeff = Wo, 



N 



. ^ 2-'= ~ V2wo. (16) 

\ k=l 



The ratio 4^ describes how well the oscillating component 



Wk exp ( J 2vklzk)lxk ew{-j2vklzk), (17) 
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is averaged. This ratio 



2vk 



2(V«Li + 1 + 1) > 4. 



(18) 



> 2 



(19) 



Vk+i Wk+i {^al+^ + 1 + 1) 



4 Non-Resonant Conditions 



We now check whether all the oscillating terms captured by Hamiltonians H'j^ and i/^" that 
were neglected result in an effective coupling. If this were to happen; in the modulation 
frame of the rf-field and chemical shifts, we will see a net coupling evolution. Therefore, we 
evaluate the evolution of the couplings in the modulation frame defined by the Hamiltonian 
Hm — Hcs + Hrf. The evolution of this frame takes the form 



such that |/n(u;)| <^ cq. Here /„(a;) represents the chemical shift in the n''^ frame. Where 



U{uj, t) = exp(-wi4j(<^)t) . . . exp{-ivkIzk{oj)t) ■ ■ ■ exp(-iv„7^„(^)i)e„(a;, t), 



Where 




(20) 



fe+i 



(21) 



where 



N 



Vk{t) = '[[exp{-ivjl,.t). 



(22) 



j=k 



In this notation, 



u{uj,t)^Vi{t)en{co,t). 
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(23) 



Where 



Qn{oo,t)=I+ f H{T)dT+ f r H{T)H{a)dTda + . . . . (24) 
Jo Jo Jo 

Then in the frame of the Hamiltonian Hm{t), the evolution takes the form 



I - i27rJ f VihS.VidT - 2nJ f [VlI,S,Vi, T H{a)dcT] dr + (25) 

Jo JO JO 

v ' 

J(r) 

where, the first integral is averaged out in the subsequent frames as — > and Vk > ^Vk+i 



in Eq. ( 19 ) prevents generation of static components from the oscillating parts. We evaluate 
the second integral in the series above. If J(r) has any static components this would reflect 
residual couplings in the system. We use the notation H'l{f) to denote the discrete set of 
frequencies present in oscillating Hamiltonian H'l{f). The discrete frequencies is the set 

{exp(-2a;/,)i/^exp(za;/,)}(/) = {uj ± Hl{f) , ±Hl{f)} (26) 



HlU) = (^kVk + E ^kVk+i. at G {±2}, hk e {0, ±1, ±2} 
Similarly, we have 

N 
k=l 

where Ck € {0,±1}. We compute the overlap of the two set of frequencies, we find the 
smallest value of |A| satisfying 

fc-l N 

j=i j=k+i 
where dj G {0, ±1, ±2, ±3}. If |A| > 0, we avoid a resonance condition. 
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4.1 Simulation and Experiments 

Simulations were carried for a carbon-proton (IS) spin system, with a coupling constant 
J = 140Hz in Eq. ([!]). The Carbon chemical shift range was Co(27r)~^ = 22.5 kHz, which 
corresponds to a 200 ppm Carbon chemical shift at 900 mHz proton frequency. At 6.25 
kHz rf-power, the | pulse corresponds to 40 fis. We choose this as our Aeff in the MODE 



sequence, so that 



3.6. We choose Wk = Wq, and calculate Wq that corresponds to 



Af,ff ~ 6.25 kHz. For = 6, this corresponds to Wq = 4.8 kHz. Fig. |2] shows the amplitude 
and phase profile of the MODE sequence. Fig. [T] shows evolution of initial magnetization 
5*^(0) = 1 under the MODE sequence for various values of N for a period of ~ 12 J^. Fig. 
[T] also shows the decoupling efficiency. 



(28) 



for various offsets. 



Decoupling PerTormance vs Modulations 



Decoupling PerTormance as Function of Chemical Shift 




time (1/J) 




Figure 1: The left figure shows how the coherence Sx{t) evolves as a function of time in 
units of J~ as the number of modulations in the Rf-field are increased. The simulations 



correspond to J = 140 Hz, A, 



eff 



6.25 kHz and Cq = 22.5 kHz corresponding to 



3.6 



and the chemical shift value = 1. The right figure shows the decoupling efficiency 



^ lo Sx{t)dt, as function of ooq. Here T 



12 

rNJ rNJ 

J 



85 ms. 
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Figure 2: The top figure shows the Amphtude (in kHz) of the rf-field as a function of time 
in units of J~, corresponding to iV = 6 modulations, with wq = 4.8 kHz. This corresponds 
to Aeff = 6.25 kHz and the maximum amphtude Amax = 12 kHz. The bottom panel shows 
the phase (in degrees) as function of time. 

The parameters of the MODE sequence for N = 6, Wq = 4.8 kHz and Cq = 22.5 kHz are 
ci = 9.1, C2 = 3.5, C3 = 1.23, C4 = .38, cg = .09 and ce = .01 kHz. The ratio |^ = 4.6, 
^ = 3.7, ^ = 2.85, ^ = 2.02, ^ = 1.26, ^ = .61 and ^ = .17. In the last frame the 
ratio of chemical shift to effective control is significantly reduced. The smallest frequency 
|A| in Eq. ([27]) is 197.3 Hz which avoids resonance. 

Fig. ([4]) shows the experimental spectra obtained on a system of methyl Iodide with 
MODE sequence. The experimental parameters of the system are the same as in simulations. 
Comparison of conventional decoupling sequences with MODE shows it is much broadband 
for same root mean square rf-power. Fig.(|5]) shows the 2D HSQC spectra of the protein 
GBl, comparing GARP and MODE decoupling sequence on ^^C during direct detection. 

The figure ^ shows the performance of the MODE decoupling sequence as function 
of rf-inhomogeneity, which is captured by the parameter e = (1 + 5), such that the actual 
amplitude of the rf field is eA where A is the nominal amplitude. 

In the absence of inhomogeneity, for all u G [— Co,Co], //c(w) < fk{co) and the ratio ^ is 
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Decouplig Perfonnance in Presence or Rrinnomogeneitv 



E = 1.1 1 - 1.2 




Q 4 I 1 1 1 1 1 1 1 1 1 

25 m 15 ID 5 5 ID 15 m 25 

Chemical shift (Khz) 

Figure 3: The above figure shows the decouphng efficiency rj, as function of Uq. Here 
T ~ y ~ 85 ms for rf-inhomogeneity values corresponding to e = .9, e = 1.1 and e = 1.2 
with 6 = —.1, 6 = .1 and 6 = .2 respectively. The mode sequence is designed with = 6, 
^ = 3.6 and cq = 22.5 kHz. 



constantly decreasing. In the presence of rf-inhomogeneity, for our simulation with 6 = .1 
and S = .2, it is observed that ^J4r^ begins to increase. For S = —.1, we observe that 
fk{uj) ^ /fc(co). The performance of the MODE sequence in presence of rf-inhomogeneity 
for 5 > is comparable to the ideal case. For 6 = —.1 a small set of offsets are affected by 



rf-inhomogeneity. For these offsets, we find that the angle nfcCos(^fc) in 12 is the largest, 
indicating residual coupling in the system. 

The following section discusses techniques to design in presence of rf-inhomogeneity. 

5 Discussion and Conclusion 

In the presence of rf-inhomogeneity, we can redefine our iterative procedure for computing the 
modulation frequencies Vk- In presence of inhomogeneity, the spread of effective shifts after 
the ffist frame transformation is from [{l — 5)wo, J {I + S^Wq + Cq], the left limit corresponds 
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to zero offset and smallest Rf-amplitude. The right limit corresponds to largest Rf-amplitude 
and chemical shift. We choose vi as center of this spread. Following this we obtain that 



'^k+i = ^ ; Cfc+i = ^ . (29) 

This gives 

ak+i = ^al + (1 + 6f - (1 - 5). (30) 

Note 



^al + {l + 5f-{l-5)>au (31) 

with equality when 

_ _ 25 

The ak values decrease until they reach Up. For 5 = .1, we have ap ~ .22, which is 
good enough as the ratio of chemical shift to Rf-strength is small. More generally, taking 
Wk+i = ^5 "we can move this fixed point close to zero. 

In this paper, we introduced the multiple rotating field technique as a means to design 
broadband heteronuclear decoupling sequences in Solution NMR. It is important to point 
out that the MODE sequences for the case when iV = 1, simply reduce to the well known 
TPPM decoupling pulse sequence pj. Also see [21 |3l HI |5|, El [7]. The main contribution 
of this paper lies in realizing the importance of the ratio and showing that by adding 
extra modulations in the rf-field and successively transforming into rotating frames, we can 
significantly reduce the ratio at and improve the decoupling performance. The technique is 
also expected to find use in designing decoupling pulse sequences in Solid State NMR and 
design of various excitation, inversion and mixing sequences, where spread of chemical shifts 
can simply be removed by transforming into a suitable K frame. 
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Figure 4: The above figures show the experimental proton spectra obtained at 500 niHz for 
a 13C-iodomethane sample, dissolved in D-chloroform (D, 99. 8%, Cambridge Isotope Labo- 
ratory, Inc.) The rf- field corresponds exactly to the simulation parameters with N — 6. The 
power level of 13.4 dB corresponds to a 40 /xs | pulse, i.e, a 6.25 kHz Rf-field. The Aeff for 
the MODE decoupling sequence is at 6.25 kHz. 
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Figure 5: The above figures show the HSQC spectra on a sample of the protein GBl obtained 
at 500 mHz . In the indirect dimension, the proton is decoupled by a refocusing vr pulse. In 
the direct dimension, the carbon is decoupled by application of GARP and MODE decoupling 
sequences with mean rf-amplitude of ~ 3.68 kHz and 6.25 khz respectively. The MODE 
sequence used is the same described earlier in the section. A plasmid pET30-a containing 
the gene coding for the HIS-tagged GBl protein (59 amino acid immunoglobulin binding 
domain from streptococcal protein G) was transformed into E. coli BL21(DE3) cells for 
protein expression. A 15N/13C uniformly labeled sample was prepared by overexpression in 
M9 minimal media containing 13C glucose and 15NH4C1. The protein was initially purified 
using a Nickel resin (Qiagen). The HIS tag was later removed with TEV protease and the 
resulting GBl was further purified on a FPLC using a Sephadex column. The sample was 
concentrated to a final concentration of 1 mM in phosphate buffer (20 mM, pH = 6.5) with 
150 mM NaCl and 100 M EDTA. 
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